Spectral Properties and Synchronization in Coupled Map Lattices 
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Spectral properties of Coupled Map Lattices are described. Conditions for the stability of spatially 
homogeneous chaotic solutions are derived using linear stability analysis. Global stability analysis 
results are also presented. The analytical results are supplemented with numerical examples. The 
quadratic map is used for the site dynamics with different coupling schemes such as global coupling, 
nearest neighbor coupling, intermediate range coupling, random coupling, small world coupling and 
scale free coupling. 
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I. INTRODUCTION 

Synchronization of large interacting systems has been 
observed in several natural situations such as synchro- 
nized flashing of the fire flies, pace maker cells of the 
heart, neurons, etc [Q, ^, [|. Synchronization of chaos 
in low dimensional systems was studied by Pecora and 
Carroll (i| . ft has also been studied in coupled oscillator 
systems and other spatially extended systems []|, [| . Due 
to potential applications in various problems of practical 
interest, synchronization of chaotic elements in a coupled 
dynamical system has been an active area of research 



Spatially extended systems are suitably modeled by 
coupled map lattices (CML). In comparison to partial dif- 
ferential equations, CMLs are more suitable for computa- 
tional studies because of the discrete nature of time and 
space while all the analytical aspects of dynamical sys- 
tems theory can also be used. CMLs were introduced as 
a simple model for spatio-temporal chaos |9) . They show 
a variety of phenomena from regular periodic behavior 
to very complicated spatio-temporal patterns, chaos, in- 
termittency, etc 10 . In CMLs, the dynamical elements 



are situated at discrete points in space, time is discrete, 
and the state variable is continuous. Each spatial unit 
is coupled to its neighbors. The selection of neighbors 
is determined by the structure of the network. In most 
studies diffusive coupling (nearest neighbor interaction) 
is used. There are studies on CMLs with various coupling 
schemes, such as open network, random network, global 
coupling , etc [pr| . In most studies a symmetric coupling 
matrix is employed. 

Here we study the synchronization properties of sys- 
tems formed by a large number of identical dynamical el- 
ements that are connected by identical symmetrical links. 
We derive general conditions for the stability of spatially 
homogeneous solutions of a CML with any symmetric in- 



teraction matrix making use of the spectral properties of 
the interaction matrix. The coupling topology can affect 
crucially the synchronizability of the system. 

In the next section we describe the properties of the 
spectrum of the CML. We perform a linear stability 
analysis and give the conditions for the stability of syn- 
chronous solutions and different regimes of stability. This 
is given in Sec. III. Results on the global stability anal- 
ysis are given in Sec. IV. In Sec. V we provide some nu- 
merical results to elucidate the analytical results with 
specific examples. Here we take a quadratic map for the 
site dynamics as an example but the results are valid for 
any dynamical system. The results are even more gen- 
eral in the sense that their validity is not restricted just 
to CMLs, but can further be applied almost directly to 
partial differential equations, coupled ordinary differen- 
tial equations, etc. Finally we provide a discussion on 
related aspects of the dynamics of CMLs. 



II. CML AND ITS SPECTRUM 



We consider a coupled map lattice of the form, 



u{x, n + 1) = e — ^ f( u (v> n )) - f( u ( x , n )) 



v 



+f(u(x,n)), 



(1) 
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where n x denotes the number of neighbors of x. Here, 
/ : R — ► R is some differentiable function, often chosen 
to be the quadratic (logistic) map in the literature, x is a 
spatial variable, its domain being some finite discrete set 
M. That set carries a neighborhood relationship, speci- 
fying which y € M are neighbors of a given x (notation: 
x ~ y). The extreme case is the one of a global coupling 
where all y are neighbors of any x. If M has the structure 
of a fc-dimensional periodic grid, the other extreme case 
is the one of nearest neighbor coupling where only those 
y are neighbors of x that are one step away from x in one 
of the coordinate directions. In that case, each x has 2k 



2 



neighbors. Of course, we also have the trivial case where 
each x is its own neighbor, but has no other neighbors. 
That case of course, represents the absence of coupling. 

In the sequel, the only assumption we shall need is that 
the neighborhood relationship is symmetric, i.e., if y is a 
neighbor of x, then x in turn is a neighbor of y. We also 
adopt the - completely inessential - convention that x is 
not considered as a neighbor of itself. (Abandoning that 
convention would simply amount to a redefinition of the 
value of e.) Finally in order to avoid trivial case distinc- 
tions, we assume that the neighborhood relationship is 
connected in the sense that for any given xi,X2 £ M, we 
find yi = xi, y 2 , . . . , y m = x 2 , s.t. y j+i is a neighbor of 
yj for j = 1,2,... , to — 1. We consider n = 0, 1, 2, . . . , 
as the time variable of the evolution. 

Our subsequent analysis will not depend in conceptual 
terms on the detailed structure of M. Of course, the 
numerical values of the bifurcation parameters below will 
reflect the geometry of M . 

Our analysis is phrased in general terms and so it is 
straightforward to extend it to the cases: 

- where / is vector valued, 

- where M is a continuous space which then has to 
carry a measure dfx, and the averaged sum needs to 
be replaced by an averaged integral, 

- to weighted neighborhoods i.e., where we are given 
a nonnegative function 

h: M x M -> R+ 

that is symmetric (h(x,y) — h(y,x) V x,y £ M) 
and consider in place of the averaged sum in equa- 
tion §) 



where \M\ stands for the number of elements of M. We 
also put \\u\\ := (it, u) 1 / 2 , (L 2 -norm of u). We consider 
the operator, 



C : L 2 (M) 



L 2 {M). 



£v(x) := — ^2 v(y) - v(x). 



(2) 



v 



C has the following properties: 

(i) C is selfadjoint w.r.t (.,.): 

(it, Cv) = (Cu, v) 

for all u,v £ L 2 (M). This follows from the symmetry of 
the neighborhood relation. 

(ii) C is nonpositive: 

{Cv,v) < 0. 

This follows from the Cauchy-Schwarz inequality, 
(iii) 



Cv = 



v = constant. 



Hence, (i) implies that the eigenvalues of £ are real. 
By (ii), they are nonpositive; we write them as —Xk, and 
the eigenvalue equation then is 

Cu k + X k u k = 0. 

We order the eigenvalues as Ao < Ai < A 2 < • • ■ < Xk- 
(This convention deviates from the one used in the lit- 
erature. Our operator C corresponds to the interaction 
matrix minus the identity matrix, and one usually con- 
siders the eigenvalues of the former in descending order.) 
We may then find an orthonormal basis of L 2 (M) , 



(the situation in (0) corresponds to the choice 

jl if x,y neighbors; 
h ^ = \o else } ' 



K) 



fe=i, 



of eigenvectors of C. 

By (iii) the smallest among the Afe is 



replacing the last term f{u{x,n)) in (|I]) by 
g(u(x,n)) for some function g, 

as well as to the case of coupled ordinary differential 
equations in place of difference equations. 



As these extensions are rather trivial, we refrain from 
carrying them out. 

The following represents a generalization of the linear 
stability analysis that has been carried out in the litera- 
ture for some special cases such as global coupling [12 , 
nearest neighbor coupling 13 , and random coupling [14|. 

We shall need the L 2 -product for functions on M: 



(u,v) 



— n x u(x)v(x), 



Ao = 0, 

and this is a simple eigenvalue (because we assume that 
the neighborhood relationship is connected), i.e., 



Afe > for k > 0. 



(3) 



lEtf 



The numerical values of the bifurcation parameters oc- 
curring below will depend only (besides on e and the 
Lyapunov exponent of /) on the eigenvalue spectrum of 
L. This eigenvalue spectrum, of course, reflects the un- 
derlying geometry of M and of the coupling. Some gen- 
eral considerations may be helpful for understanding this 
point. 

In the case of global coupling (including self coupling) , 
we have 
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Ao = (as always) 



and 



since 



Ai — A2 



A fc = 1, 



Cv = —v 

for any v that is orthogonal to the constant map, i.e., 
satisfies 



y&M 

If we shrink the neighborhood size, then the eigenvalues 
can separate and grow, and in particular, the largest one, 
Xk, will become larger the smaller the neighborhood size 
is. In particular, 

Xk > 1 

as there may exist v 6 L 2 (M) with 

xeM v 

(e.g. M — {1, 2, . . . , m}, m even, m > 2, with /1 having 
neighbors fi—1 and closed periodically, i.e. m+1 = 
1, 



1. 



even, 
odd 



Conversely, if the neighborhood interaction matrix of 
all points is the same and kept fixed while we increase 
the size of M, then all eigenvalues will decrease. This is 
a version of Courant's monotonicity theorem Jl5| |. Thus, 
from our analysis below, synchronization will require, if 
possible at all, a larger value of the coupling parameter 
e. 

We also have the following version of Courant's nodal 
domain theorem [jisfl : 

Lemma 1. Consider M as a graph Tm, with an edge 
between x and y precisely if x and y are neighbors. Let 
Uk be an eigenfunction for the eigenvalue Xk, with our 
above ordering, = Ao < Ai < A2 < • • • < Xk- Delete 
from the graph Tm all edges that connect points on which 
the values ofuk have opposite signs. This divides Ym into 
connected components T±, . . . , T/. Then I < k + 1. 



III. LINEAR STABILITY ANALYSIS 

We now consider a solution u(n) of the uncoupled 
equation, 



u(n + 1) = /(«(n)). 



(4) 



Clearly, u(x,n) = u{n) then is a solution of (|). This 
solution is spatially homogeneous, or as one says, syn- 
chronized. The synchronization question then is whether 
for certain values of the coupling parameter e, any so- 
lution of (Q) asymptotically approaches a synchronized 
one. A somewhat weaker question is whether, when we 
consider a perturbation 



u{x,n) = u{n) + 6ctk(n)uk{x), 



(5) 



by an eigenmode Uk for some k > 1, and small enough 
S, afe(n) goes to for n — ► 00, if u(x, n) solves (0). That 
question can be investigated by linear stability analysis 
and we proceed to carry that out. Inserting (H) into (|l|) 
and expanding about 5 = yields 



Q.k{n + !) = (!- eX k )f'(u(n))a k (n), 



(6) 



/' denoting the derivative of /. So the sufficient local 
stability condition 



hm — log — r- 

n^oo n 6 a k (Q) 



becomes 



JV-l 



lim log TT 



ctkjn + 1) 
a k (n) 



< 



(7) 



N-l 



log|l-eA fe |+ lim l^log|/'(u(n))|<0. (8) 

N—>oc 1\ L — ' 



n=0 



Here, 



N-l 



N-><x N 

n=U 

is the Lyapunov exponent of / and so the stability con- 
dition (ph is 



We may have 



|e"»(l-cA fc )| < 1. 



/i > 0, 



(9) 



(10) 



i.e. temporal instability, but (^|) for all k > 1; i.e. syn- 
chronization. We shall now assume ([Io|) for the remain- 
der of this section. By our ordering convention for the 
eigenvalues, (H) holds for all k > 1 if 



1 - e-^° 1 + e-^° 

< e < 



Ai Xk 
In order to satisfy that condition, we need 

Xk_ d l ° + 1 
Ai e^o - 1 ' 



(11) 



(12) 



By our above discussion this hold in the globally coupled 
case because there Xk = Ai. By way of contrast if we 
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have nearest neighbor coupling, this can only hold if the 
size of M is not too large. (For a 1-dimcnsional chain, 
the critical size is 5, with a large value of e. If we have 
second nearest neighbor coupling, the critical size of a 
one dimensional chain is 9.) 

Let us now assume that (|l^) holds. We then pre- 
dict the following behavior of the coupled system as e 
increases. 

For very small values of e > 0, as we assume ( J]~0| ) 

e"»(l-eA fc ) > 1, 

and so, all spatial modes itfc, k > 1, are unstable, and no 
synchronization occurs. If we are in the globally coupled 
case, then there exists a single critical value e c such that 

e"»(l-e c A fc ) = 1 

for all k = 1, 2, . . . , K, For e > e c , the dynamics become 
synchronized. For e slightly smaller than e c , one observes 
intermittent behavior, clustering, etc |l7| . 

Let us now consider the more interesting case where 
the coupling is not global so that not all the X k are equal; 
in particular 

Ai < X K . 

We then let e& be the solution of 

e^(l~e k X k ) = 1 

The smallest among these values is 6k, the largest e%. If 
now, for k\ <k 2l 

€k 2 < e < e kl 

then the modes Uk 2 ,Uk 2 +x, ■ ■ ■ ,Uk are stable, while the 
modes tti,lt2, . . . , u kl are unstable. Because of Lemma 1, 
we see that desynchronization can lead to utmost k 2 + 1 
subdomains on which the dynamics is either advanced or 
retarded. 

In particular, if e increases, first the highest modes, i.e., 
the ones with most spatial oscillations, become stabilized, 
and the mode u\ becomes stabilized the last. So if e 2 < 
e < ei, then any desynchronized state consists of two 
subdomains. 

We then let e k be the solution of 

e^(e- fe A fe -l) = l 

Again, 

Because of (p"l|), 

ei < l K - 

If 

ei < e < e K , 



then all modes u k , k — 1,2,... , K, are stable, and the 
dynamics synchronizes. 

If e increases beyond ex, then the highest frequency 
mode uk becomes unstable and we predict spatial oscil- 
lations of high frequency of a solution of the dynamics . 
If e increases further then more and more spatial modes 
become destabilized. 



IV. GLOBAL STABILITY ANALYSIS 

The basis of the preceding analysis was a linear ex- 
pansion about a synchronized state u(n). Therefore, that 
analysis is valid only for small perturbations about such 
a state. In this section, we want to derive a criterion that 
guarantees synchronization for arbitrary starting values 
m(x,0) of a solution of (Q). 

From general principles of functional analysis (see 
fill), there exists an operator, 

A : L 2 (M) -» L 2 (M) 

with 

~{u,Cv) = (Au, Au),Vu,w e L 2 (M). (13) 

This follows from the self adjointness of C. It is not diffi- 
cult to write a A down explicitly, but our more abstract 
approach provides the advantage of a less cumbersome 
notation. 

A is nonnegative in the sense that 

(Am, Aw) > 0,Vu G L 2 (M), (14) 

and we even have 

Au = u = constant (15) 

(This follows from the nonpositivity properties of C). 
Moreover A commutes with C, i.e., 

AC = £A, (16) 

and so, we may assume that the u k are also eigenfunc- 
tions of A. 

Therefore a natural ansatz for a Lyapunov function for 
the dynamics ([!]) is 

$(n) := (Am(., n), Au(.,n)), (17) 

and it remains to derive conditions under which 

$(n) — » 0, for n — > 00. (18) 

We have 

$(n + l) = (Au(.,n + 1), Au(.,n + 1)) 

= (A«(., n + 1), A(eL/(., n) + (1 - e)/(., n))) 

by ©. 
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Since the u k are an orthogonal basis of L 2 (M) , we may 
write 



In the special case of global coupling, the synchroniza- 
tion condition becomes, 



A" 



f(u(x,n)) = ^2l3 k (n)u k (a 



k=0 



with /3fc(n) = (/ (u(. , n)) , u k ) ■ Inserting this into the last 
equality, we get 

$(n+l) = (a U (., n + 1), A^(l - e\ k )j3 k (n)Uk \ . 

V k=0 / 

(19) 

The important observation now is that in the last sum, we 
can discard the summand k = 0, because uo is constant, 
and so 

Aw = 0. 

Moreover, we observed above that, since A commutes 
with £, we may assume, 



(Ait*,, Aui) = 0,for k ^ I, 



and so 



||A/(u(.,n))|| 2 =^^(n)||A^|| 2 . 

k=Q 

Using these observations and the Cauchy-Schwarz in- 
equality in Jis|), we may estimate 

$(n + 1) < ~||Au(.,n + 1)|| 2 + 5(1 - eA a ) 2 ||A/( W (.,n))|| 2 , 

(20) 



assuming |1 — eA^-l < 1 — eAi, i.e., 

2 

e < 



(21) 



Ai + A A - 

If we now use the coarse estimate 

||A/(u(.,n))||<Bup|/'|||Au(.,n)||, (22) 
we obtain from (|20| ) 

$(n + 1) < (1 - eAi) 2 sup |/'| 2 $(n). (23) 
We conclude 

Theorem 1. The coupled dynamical system asymp- 
totically synchronizes if e satisfies (pJp and 



(l-eA 1 )sup|/'| < 1. 



(24) 



Remark: If ( pl| ) does not hold, (|24| ) needs to be re- 
placed by, 



(l-e)sup|/'| < l,(0<e< 1) 



(26) 



(eA K -l)sup|/'| < 1. 



(25) 



The reason why we have sup \ f'\ in (pi} ) , in place of 
e^°,/io being the Lyapunov exponent of /, as in section 
3, is that here we do not linearize about a spatially ho- 
mogeneous solution. Our global approach rather requires 
to consider any solution u(x, n) of (|l|). This means, how- 
ever, that our condition ([24]), while sufficient, need not 
be necessary for synchronization. 



V. NUMERICAL RESULTS 

In this section we demonstrate our results with dif- 
ferent coupling schemes or network topology. For our 
numerical study we took the quadratic map for the site 
dynamics. The quadratic map is a widely studied chaotic 
map, given by f(x) = 1 — ax 2 fisfl . Here a is a param- 
eter and varying its value the single map shows a vari- 
ety of dynamical phenomena. It becomes chaotic when 
a « 1.4011, going through a period doubling bifurcation 
sequence. At a = 2 the map is maximally chaotic, with 
a Lyapunov exponent /j,q — log(2). 



A. Global coupling 

In the case of global coupling, we have Ao = and 
Ai = A2 = ■ • • = A m _i = 1. (The self coupling term 
is also included here.) This case has been studied in 
various contexts. When exp(/xo)(l — e) < 1, the spa- 
tially homogeneous solution is stable, as shown in [fl2[ . 
For the quadratic map with a = 2, it becomes sta- 
ble when e > 0.5. Just below this value the system 
shows spatio temporal intermittency, clustering phenom- 
ena, etc |l(}. In Fig. |l|(a), we display er(n), the fluc- 
tuation of the state variable from the mean, defined by 
a 2 (n) — ^^2™= 1 {xi(n)—x(n)) 2 ,(x(n) is the average of all 
Xi(n)), for different values of e, for the case with a = 2.0. 
It can be seen that when e > e c = 0.5, the value of a 
becomes zero (within the numerical accuracy) indicating 
that the system is synchronized. Though the linear sta- 
bility does not guarantee the synchronization from ar- 
bitrary initial conditions, in this case it happens. We 
started with random initial conditions for the individual 
sites, and after a few iterations the system synchronizes, 
indicating the stability of the spatially homogeneous solu- 
tions in these parameter regimes. For a — 1.9, the system 
synchronizes for a smaller value of e, since the Lyapunov 
exponent at that parameter value is 0.5490. Here the 
critical value is e c = 0.4225. Figure |l|(b) gives details of 
this case. We took m — 1000 for our simulations. 
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FIG. 1: Fluctuation of the mean, <r(n) is shown as a function 
of the coupling strength, e. At each value of e, 200 final 
iterates of o(n) axe plotted. Here m — 1000 and the coupling 
is global. In (a) a = 2.0, and in (b) a = 1.9 . 



FIG. 2: Same as in Fig. |l| with nearest neighbor coupling, for 
m = 5 , in (a) a = 2.0, and in (b) a — 1.9. 



B. Nearest Neighbor Coupling 

Here the eigenvalues are given by Ao = and A; = 
1 — cos( — ) i = 1,2,3,... ,m — 1. The first nonzero 

V 771 ' ' 111 1 

eigenvalue is 

Ai = 1 — cos( — ), 
m 

and the largest eigenvalue is 

I 2 for even m 

|l + cos(— ) for odd m 

Using this one can calculate the maximum value of 
m at which the spatially homogeneous solution can be 
stable using the condition for linear stability. It will occur 
when 

A_g cxp(^p) + 1 
Ai exp(^o) - 1 
and the value of e lying between 

1 - exp(-/j ) <e< 1 + exp(-/x ) 
Ai \k 



For a CML with a fully chaotic quadratic map the 
maximum value of the system size which can sustain a 
stable synchronous solution is m = 5, when e is between 
0.7236 and 0.8292. In the case of m = 6, the first mode 
becomes stable at e — 1, but the last mode becomes 
unstable for a value of e above 0.75. Hence there is no 
synchronization. The second mode is stable when e is 
between 0.333 and 1. In Figs. we give the plot of 
the fluctuation of the mean field for different values of e, 
for rn = 5 and m = 6, when a = 2.0 (a) and a = 1.9 
(b). When m = 6, between e = 0.33 and 0.75 only one 
mode is unstable. From the spectrum one can see that 
the largest value of m for which only the first mode is 
unstable is m = 9 when 0.72 < e < 0.75. For higher 
values of m more than one mode will be unstable for any 
value of e. So there cannot be synchronization in large 
systems with nearest neighbor coupling. 



C. Intermediate Range Coupling 

If we consider k nearest neighbors (there will be 2k 
neighbors for each site) the eigenvalues are given by Xq — 
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FIG. 3: Same as in Fig. |^ with nearest neighbor coupling, for 
m — 6 . 



FIG. 4: Same as in Fig. [| with two nearest neighbors coupling, 
for m = 9. 



and 



1 y ^ 2irij 
X * = 1 -^l^ C0S (^r)' * = 1, 2, . • • , m - 1. 



Let us consider the case of two nearest neighbors 
(k = 2). As in the case of NN coupling one can find 
the maximum value of m at which the CML can sustain 
stable synchronous chaotic oscillations. For k=2, it is 
m = 9 with 0.33 < e < 1. The maximum value of m at 
which the second largest mode also becomes unstable is 
m = 18. Figure gives the plot for m = 9 for a = 2.0 
(a) and a = 1.9 (b). For 3 nearest neighbors (k = 3), it 
is at m = 12, and for fc = 4, m = 15. 

One can see that for k/m > 0.301 the system syn- 
chronizes when the coupling is strong, i.e., e = 1. In 
Fig. ||, cr(n) is shown as a function of k for m = 1000 
near the synchronization transition region. In each grid 
corresponding to a k value, 1000 final iterates of <r(n) 
are plotted after discarding initial transients. The sys- 
tem shows synchronization when k = 301. As the system 
size increases we need a higher number of neighbors for 
synchronization. For a fixed number of neighbors the 
behavior is like that of the NN case; there is no synchro- 
nization when the system size increases. 
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FIG. 5: Here a(n) is plotted for different values of k (denoted 
inside the grid) with intermediate range coupling, for m = 
1000, a = 2.0, and e = 1.0. Between two grid lines 1000 
iterates of cr(n) are plotted. 



D. Random Coupling 

Now we consider a case where there is coupling between 
random sites. For every site we randomly select k other 
distinct sites and connect them with each other under the 
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FIG. 6: Same as in Fig. g with random coupling for different 
k values, for m — 1000, and e — 1.0. In (a) a = 2.0 and in b) 
a = 1.9. 

constraint that self and multiple coupling is prohibited. 
The average degree of a node in such a graph obtained is 
2k. For the quadratic map with a = 2.0 and e = 1, the 
system synchronizes for large m, if k > 8, in contrast to 
the unsymmetric case where it does so for k > 4 [1^] . 

We plot the fluctuation of the mean field, ex, for differ- 
ent values of k for to = 1000, e = 1.0, a = 2.0 (Fig. |(a)), 
and a = 1.9 (Fig. ||(b)). It can be seen that the system 
synchronizes when the average degree of a vertex is 8 or 
more, for the completely chaotic quadratic map. This is 
independent of the system size to. From random matrix 
theory one can see that the value of Ai depends only on 
k f20| . For smaller to, synchronization can occur below 
k = 8 because of the finite system size effects. So unlike 
in the case of nearest neighbor or intermediate range in- 
teractions, in the case of random coupling, one can have 
chaotic synchronization for any arbitrarily large value of 
to, if the number of neighbors (k) is larger than some 
threshold determined by the value of the maximal Lya- 
punov exponent of the chaotic map. 



E. Small-world Networks 

Small-world (SW) networks have an intermediate con- 
nectivity between regular and random networks. They 
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FIG. 7: Same as in Fig. |l|with small- world coupling, for m = 
1000, k = 10, a = 2.0, and p = 0.1. 
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FIG. 8: u(n) for different values of p with small-world cou- 
pling, for m = 1000, k = 10, a = 2.0 and e — 1. Between two 
grid lines, corresponding to a p value, 1000 iterates of cr(n) 
are plotted and the corresponding log 10 (p) is denoted at the 
bottom. 



are characterized by a very small mean path length as in 
random networks while at the same time having a high 
clustering coefficient as in regular networks. SW coupling 
is done as in the Watts and Strogatz algorithm ]2l]| . We 
start with a lattice of m vertices each connected to its k 
neighbors. With a probability p we reconnect each edge 
to a vertex chosen uniformly at random over the entire 
lattice. Duplicate edges are avoided. It has been shown 
in pT) that even for a very small random rewiring prob- 
ability p there is a transition to the small-world regime. 

Here we took p = 0.1 at which there is small world 
effect on the structural properties of the graph. Figure 
gives the fluctuation <j(n) for different values of e, k = 10 
and to = 1000. One can see that there is no synchroniza- 
tion at this value of p, \i « 0.08 . When p = 0.8, there is 
synchronization for e = 1.0. At this value the number of 
random connections per vertex reaches the value needed 
for synchronization. In Fig. H, u(n) is plotted for p val- 
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FIG. 9: Same as in Fig. |lj with scale free coupling, for m 
1000, k = 6 and a = 1.9. 



ues from 0.0001 to 1, and e = 1.0. Between two grid lines 
1000 iterates of a(n) are plotted and the corresponding 
log 10 (p) values are also denoted. From this figure we can 
easily see that there is no synchronization for smaller p 
values. 



and to be independent of the time of evolution (size) of 
the network. 



VI. CONCLUSION 

We studied the spectrum of coupled map lattices and 
its relation to the stability properties of the spatially ho- 
mogeneous solutions. We derived conditions for the ex- 
istence of such solutions using linear stability analysis. 
Conditions obtained from a global stability analysis are 
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FIG. 10: Same as in Fig. |^ with scale free coupling, for 
1000, a = 1.9 and e = 1.0. 



Another widely studied class of networks are the scale 
free networks, where the degree distribution obeys a 
power law which is observed in many real networks. We 
studied the synchronization of a scale free network con- 
structed by the Barabasi- Albert algorithm Q . We start 
with ko vertices and at every time a new node is intro- 
duced. The new node is connected to k already existing 
nodes and they are selected with a probability propor- 
tional to the degree of that node. The process is contin- 
ued for a long time and then the degree distribution is 
described by the power law, P(k) ~ k~ 1 , where 7 = 3. 
It is independent of kg. For this study we took fco = k 
and a network of size m — 1000. Figure ^ shows the a 
versus e plot for k — 6, and a — 1.9. In Fig. [l0[ a(n) is 
plotted for different values of k, for the case a = 1.9. The 
synchronization behavior is comparable to that of a ran- 
dom network. When k > 8 there is synchronization for 
a = 2.0. We checked our results with higher values of m 
also. The results seem to converge for large system sizes 



also provided. Our results are supplemented with nu- 
merical examples. For the numerical study the quadratic 
(logistic) map is used for the site dynamics. We studied 
the synchronization properties of coupled map lattices 
with different coupling topologies such as global coupling, 
nearest neighbor coupling, intermediate range coupling, 
random coupling, small-world coupling and real-world 
coupling. The coupling topology can crucially influence 
the synchronizability of the CML. Our study can be gen- 
eralized almost directly to other spatially extended sys- 
tems. 
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